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مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری ترمومکانیکی و مکانیکی * 
علی جمالو "" . . مهدی خداپرست ۲ سیدعلی قریشیان امیری"" 


چکیده در سال‌های انعیر مدل‌های رفتاری متعددی برای پیثر‌بینی رفتار نعاک تحت اثر بارهای ترمومکانیکی ارائه شده‌است. از سوی دیگر 
برر سی رفتار شمع‌های حرارتی و حاک مجاو رآن تحت بارهای ترمومکانیکی و باا ستفاده از مدل‌های رفتاری که د رآنها تغییر رفتار حاک در 
ثر تغییرات دما درنظ رگرفته می‌شود, از اهمیت زیادی برنحوردار است. در این تحفیق یک مدل رفتاری ترمومکانیکی به نرم‌افزار المان محدود 
آباکوس | ضافه شدها ست تا ا زآن برای برر سی عملکرد شمع‌های حرارتی, | ستفاده شود. نتایج بهد ستآمده از مدل سازی عدد یآزمایش 
سه‌محوری نشان می‌دهد که کد نوشته‌شده برای مدل رفتاری مورد نظر توانایی پیش‌بینی رفتار حاک در شرایط ترمومکانیکی را دارد. هم‌چنین 
نتایج به‌د ستآمده از تحلیل شمع حرارتی باا ستفاده از مدل‌های رفتاری مکانیکی دراک ر-پراگر, کم‌کلی | صلاح شده و ترمومکانیکی تورچی 
بیانگ رآن است که باوجود تأثیر قابل‌توجه انتحاب مدل رفتاری بر نایج در تحلیل شمع حرارتی» تفاوتی میان انتخحاب مدل رفتاری مکانیکی و 
ترمومکانیکی برای نحاک کنا رآن وجود ندارد. 


واژه‌های کلیدی شمع حرارتی. مدل رفتاری خاک. ترمومکانیک اثر دما بر خاک آباکوس. 


0ص آم‌اصمل۱۷۲6-مصصمص اه وا زرا معصهصر۳۵)0 ممازظ ومهص اه ما اتمه ۱ 16 
۷۵6 1600۲6) 06 [2ع1صفملعع]۱۷ 


0 .۸ تتصم صهتطوتهمدت ۸ ٩.‏ 0202۵ 1۷۲۰ 


۲ ۵۵1۱0۷۲0۳۴ ]501 ۲۵۵/61 10 0۳0۵05۵0 0۵۵۱ ۵۷۵ ۱۱۱۵015 60۱ 6۷۵۵ ۵۵۲۵ ۳۵۵۵۵ ۱ ]۸5۳9۵6 
۳ 07 ار ۵۱۵۲۵ 0۴ ۵/۵۷۵ ۱6 رای 1۵ ۱۱۱۳۵۵۳۱۵۲۱ 15 31 ,۱۵۱۱۵ 01/۱6۳ 6 6۵۱۱ .۱۵۵05 171611۱۱0۳۱۵6۱۵۵ 
٩0۵۱ ۲‏ ۱۷۵/۱۲۵۱ 1 ۱۱۵۵615 60۳۵۹۳۵۲۷۵ 98۵ را ۵۵ هه آمع۱۱۵۳۱۱۵۱۱۵۵۵ ۱۱۵ کازهی 00/060 
۱۱۵۵۵ 60۱۱۷۵ 1۱۵۳۱۱۱۵۱۱۵۵۱۵۵۵ ۵ راک کز:ا ۱ 60۳۵۵۳۵4 0۳۵ 0۱0 1۵۴۵۵۳۵۵۲۵ 0610 6008۵65 
۰ 616۲8 0۴ 6101111۵۳۱6۵ 1116 6۷۵/۱۵1۵ 10 560 ۵ 10 5۵/۷۵۲۵ 6161161 1116 ۸۵5 1۱2 10 00060 ۱۷۵5 
1 1116 ۲۵۲ ۷۷۲۲۲۵۲۷ 6006 1۱6 ۱۵۲ مدای 1651 ]1۳۱۵2/۵ ۵ ۱۱۱۵۵61۱۵ ۱۱۵۲۱۵۵ ۲۵۷۱ 01۵1۱۱۵0 ۳۲65۱06 
6 ,۸50 .60۳0۱0۱5 ۵۳۱۱۵۱۱۵۵۱۵۱6۵ ۱۵6۳ ۵09۵۲ آز50 2۲۵۵۱6۱۱8۵ 0۴ 0۵602 کذ ۱۵۵1 60115111۷ 
,01۳۱۵۵۰۳۵۵۵۲ داع۱۱۱۵04 60۱0۲6 ۱۱6۵۱۵۱۱۱6۵ ها ار 2168 0۴ داع06 16 ۵۸۲ 0۳21۵۱۱۱۵۵ ۲۵۱۵15 
ره ۱/99 
016۲۵۷6۵ ۱۵0 5 ۱۱۵۳۵ 9 ,۲6۵۵6 ۱6 0۴ ۱۱۵۵61 0و5 ۵ دم ۱۱۵ 60۳/۱۷۲۵ ]0 6110166 1116 
011 0006610 ۱۵ ۵۲ ۱۱۵۵ 60۵ ۱۱۵۱۱۱۵۱۱۵۵۱۵ 0 ۱۱۵۵۱۵۱۵۵ ۵ 6005۵ 


٩011, ۰‏ 08 ]61160 متااجتهمره ] رمتصهطهم مور مصمط ز رلزمع ۵۶ [م0مظ متانتافومن بماژم ۳۲۵۵۵1 ۲۷۷۵۳0 م1 


* تاریخ دریافت مقاله ۹۸/۷/۷ و تاریخ پذیرش آن ۹۹/۸/۳ می‌باشد. 
(۱) دانشجوی کارشناسی ارشد. گروه مهندسی عمران دانشگاه قم. 
(۲) نویسنده مسئول: دانشیار گروه مهندسی عمران دانشگاه قم. 0061۳ 1002021256۵ :211 


۳ پژوهشگر ارشد گروه مهندسی ژئوتکنیک دانشکده مهندسی, دانشگاه علم و صنعت نروژ. 


مقد مه 
در سال‌های اخیر در مناطق مختلف حهان تیان 
ساختمان‌های مسکونی و تجاری با رویکرد بهینه سازی 
م صرف انرژی باا ستفاده از شمع‌های حرارتی تبدیل به 
یک فناوری نوظهور شده‌است [1] لذا توجه به 
تحقیقات حوزه‌ی ژئوتکنیک در این خصوص که به 
بررسی تغییرات مقاومت و تغییر شکل خاک ناشی از 
تغییرات دما پرداختهاست. بیش از گذشته مورد توجه 
قرار گرفته‌است [2]. 
مطالعات گذشته در این زمینه به دو دسته مطالعات رفتار 
خاک دحت اثر دما و م طااعات عملکرد شمع های 


حرارتی قابل تةسیم است. برا ساس تحقیقات انجام شده 


فوتای ول افرانتن ماس اک فص به کاهی فش 
پیش تحکیمی و شیب خط حالت بحرانی خاک می شود 
[3,4]. هم‌چنین محققان دیگری اثر دما بر مقاومت پرشی 
را بررسی کرده و نتایج متفاوت و متضادی از جمله 
فزایش مقاوفت برشتی [9-7] و کاهش بقاوفت فزشتی 
[5,8,9] را گزارش کرده‌اند. هم‌چنین گزارش‌های برخی 
ز مطالعات گذ شته تأثیر ناچیز دما بر مقاومت بر شی را 
نشان داده‌است [6,10], 


ز طرفی در سه دهه‌ی اخیر مدل‌های رفتاری 
ترموم کانیکی مت عددی برای خاک های مختلف ارا ثه 
شده‌است که برای رس‌هاء مدل‌های رفتاری برمبنای 
کم‌کلی اصلاح‌شده سهم قابل توجهی از آنها را 
دربر گرفته‌است [4,11-15]. عملکرد چند مدل رفتاری 
ارائه‌شده برای خاک‌های رس اشباع تا سال ۲۰۱۳ توسط 
هانگ و همکاران [16] مورد بررسی قرار گرفته‌است. 
برر سی‌های وی ذشان‌دهنده این است که تمامی مدل‌ها 
پیش‌بینی معقولی از رفتار ترمومکانیکی خاک رس اشباع 
داشته‌اند. با این حال هر مدل رفتاری دارای محدودیت‌ها 
و فرضیاتی در تئوری خود هستند که در نتایج موثر 
بوده‌است. هانگ و همکاران [17] مدل ارائه‌شده توسط 
کوی و همکاران [4] را گسترش دادند و دو مکانیسم 
پلاستیک 1۷ و را به‌ترتیب برای انقباض و انبساط 


نشريه مهندسی عمران فردوسی 


مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


حرارتی در نسبت پیش تحکیمی بالا و اثر دما بر فشار 
پیش تحکیمی معرفی کردند. درنهایت به‌منظور بررسی 
عملکرد مدل رفتاری ارائه‌شده نتایج آزمایشگاهی یک 
آزمایش سه‌محوری زهکشی شده برروی خاک رس بوم 
با نتایج به‌دست‌آمده از مدل رفتاری مقایسه شد. 

شیونگ و همکاران [2] نیز براساس مدل رفتاری کم‌کلی 
اصلاح‌شده یک مدل رفتاری ترمومکانیکی برای خاک 
ارائه داده‌اند که با قرار دادن مقدار 0607۴ برابر یک در 
آن, مدل با مدل رفتاری کم‌کلی اصلاح‌شده برابر می‌شود. 
طبق بررسی‌های انجام‌شده توسط این محققان رفتار 
ترمومکانیکی خاک وابسته به تأثیر تغیبرات دما بر دو 
پارامتر تتش پیش تحکیمی و نسبت تنش بحرانی است و 
بررسی تعادل بین آثر تغییرات دما بر این دو پارامتر 
امکان تعیین رفتار سخت شونده يا نرم شونده‌ی حرارتی 
را ممکن می‌کند. هم‌چنین مطالعات گسترده‌ای درزمینه‌ی 
برر سی عملکرد شمع‌های حرارتی به صورت تک -18] 
[21 و گروه [1,22-24] انجام شسده‌اسست. لالویی و 
همکاران [۱۸] در یک مطالعه‌ی عددی و صحرایی به 
بررسی عملکرد شمع‌های حرارتی پرداخته‌اند که در 
مدل‌سازی عددی آنها از یک مدل رفتاری 
ترموالاستوپلاستیک برمبنای مدل رفتاری دراکر-پراگر 
استفاده شدها ست. طبق تحقیقات ای شان تنش محوری 
ناشی از بارگذاری حرارتی و بارگذاری ترمومکانیکی 
بزرگ‌تر از تتش محوری نا شی از بار مکانیکی می‌با شد 
و افزایش دما به‌میزان یک درجه‌ی سانتی‌گراد ۱۰۰ 
کبلونیوتن بار محوری را افزایش می‌دهد. 

آماتیا و همکاران [25] مکانیسم ارائه‌شده را در دو حالت 
گیرداری کامل و آزاد بررسی کردند. نتایج آنها نشان 
می‌دهد که تنش نا شی از بار حرارتی درحالتی که دو سر 
شمع به‌طور کامل گیردار باشد از سایر حالات بیشتر 
هه تا هد اش یشان فان دهتنه افت ای تسا ون 
بارگذاری ترمومکانیکی نسبت به بارگذاری مکانیکی نیز 
شابن آکروچ و همکاران [26] یک آزمایش کشش 
شمع برروی شمع حرارتی در خاک با پلاستیسیته‌ی بالا 


سال سی و سوم شمارة سهء ۱۳۹۹ 


مهدی خحداپرست- علی جمالو - سیدعلی قریشیان امیری 


انجام داد ند. نتایج آز مایش های آذها بیان می کند که 
افزایش دمای خاک منجر به افزايش نرخ خزش می شود 
به‌طوری که باتوجه منحنی بار- جابه‌جایی» نشست شمع 
حرارتی در طولانی مدت (۵۰ سال) حدود ۲,۳۵ برابر 
نشست شمع عادی در همین شرایط می‌باشد. 
خدایرست و همکاران [20] به مطالعه‌ی عددی اثر دما 
بر ظرفیت باربری و نشست تحکیمی شمع‌های حرارتی 
در خاک‌های ریزدانه‌ی اشباع پرداختند و مدل رفتاری 
دراکر- پراگر را برای پیش‌بینی منحنی تنش و کرنش 
خاک درنظر گرفتند. تحقیقات ایشان نشان داد که 
افزایش مدول الاستیسیته. نسبت پواسون و نسبت 
تخلخل تحت بار گذاری حرارتی مو جب میزان کاهش 
بیشتر در تنش موثر قائم درمقایسه با شرایط ژئواستاتیکی 
مق تبود: 

نگویان و همکاران [27] به بررسی رفتار باند مدت 
شمع های حرارتی بااستفاده از یک مدل فیزیکی 
پرداختند. در تحقیق آنها بارهای مکانیکی مختلفی (از ۰ 
تا ۲۰ درصد مقاومت شمع) بر آن وارد شد و هم‌زمان با 
آن ۰ سیکل حرارتی نیز به شسمع اعمال گرد ید. 
برر سی‌های آنها ذشان داد که اعمال سیکل‌های حرارتی 


و اعمال بار ثابت منتهی به نٌ ت پلاسستیک 


(برگشت‌ناپذیر) در سر شمع می‌شود. تورچی [28] یک 
مدل رفتاری ترمومکانیکی برا ساس مدل رفتاری کم‌کلی 
اصلاح‌شده و نظریه‌ی حالت بحرانی. برای خاک رس 
در شرایط اشباع ارائه داده‌است که در آن کرنش‌های 
الانستیک ناشسی از افوایشن دماه اثر دما بر نیب ال 
بحرانی و تأثیر نسبت پیش‌تحکیمی در مدول برشی. 
درنظر گرفته شدها ست. پارامترهای | ضافه شده به مدل 
کم‌کلی | صلاح شده در تحقیق آنها باا ستفاده از آزمایش 
سه محوری به‌آسانی قابل دست‌یابی بوده‌است. آنها پرای 
صحت سنجی» نتایج مدل رفتاری خود را با ۵ سری 
داده‌های آزمایش‌گاهی موجود در ادبیات فنی کی تا 
مقاٍ سه کرده‌اند که توانایی منا سب مدل رفتاری آنها در 


پیش‌بینی رفتار ترمومکانیکی خاک های ر سی | شباع در 


سال سی و سوم شمارة سهء ۱۳۹۹ 


۵۱ 


این مقایسه قابل مشاهده‌است. جزئیات فرمولاسیون مدل 
رفتاری آنها در بخش بعد به‌طور کامل آمده‌است. 

در این مقاله عملکرد شمع حرارتی در خاک رس 
اشباع تحت بارگذاری حرارتی و مکانیکی بااستفاده از 
سه مدل رفتاری متفاوت بررسی شده‌است. در مطالعات 
عددی کل که عم نا از مدل‌های رفتاری مکانیکی و 
ترمومکانیکی متفاوتی برای توصیف رفتار خاک مجاور 
شمع حرارتی استفاده شده‌است. بر این اساس هدف از 
مطالعه‌ی حا ضر تعیین میزان تأثیر مدل رفتاری خاک بر 
نتایج به‌دست‌آمده از مدل‌سازی شمع‌های حرارتی 
بوده‌است. برای این منظور ابتدا یک سابروتین یومت 
برای مدل رفتاری تورچی [28] که برای پیش‌بینی رفتار 
ترمومکانیکی خاک رس | شباع ارائه شدها ست نو شته 
شد. سپس برای صحت سنجی کد مذکور از یک مدل 
سه‌بعدی آزمایش سه‌محوری تحکیم‌یافته‌ی زهکشی‌شده 
0) استفاده گردید و نتایج تست | لاه با نتایج 
آزمایش‌گاهی موجود در ادبیات فنی مقایسه شد. 
درئهایت نیز یک مدل سه‌بعدی المان محدود برای 
تحلیل عملکرد شمع حرارتی بااستفاده از مدل رفتاری 
تورچی [28] و مدل رف تاری دراکر- پراگر و کم‌کلی 
اصلاح‌شده برای شبیه‌سازی رفتار خاک در آباکوس 
۷ ساخته شد. هم‌چنین نتایج به‌دست‌آمده از این دو 
یتنس خی را مشاه ماش شاه عزو 
لالویی و همکاران [18] مقایسه گردیده‌است. 


مدل رفتاری و مدل‌سازی عددی 
مدل رفتاری 
در این تحقیق برای توصیف رفتار ترمومکانیکی خاک از 
سابروتین پومت مدل رفتاری تورچجی [28] که در طی 
این مطالعه نوشته‌شد. استفاده گردیده‌است. این مدل 
رفتاری با گسترش مدل حمیدی و خزایی [15] به یک 
مدل ترموالاستوپلاستیک تبدیل شلده‌است و در آن 


ایده‌ی تفاوت میان منحنی تحکیم ایزوتروب در دمای 


نشریه مهندسی عمران فردوسی 


۲ 


مرجع و دمای بالا می‌باشد. هم‌چنین در مدل تورچی 
[28] تأثیر تغییرات دما بر شیب خحط حالت بحرانی, 
سطح تسلیم و نسبت پیش تحکیمی اضافه شده‌است که 
در مدل حمیدی و خزایی [15] درنظر گرفته نشده‌بود. 

روابط مدل ارائه شده در فضای تنش سه‌محوری از 
دو مت ال کش باتک تفگ نکن ک وی 
قسمت الا ستیک برای محا سبه‌ی کرنش‌های حجمی از 
زانط و و( فاد شلات 


)۱( وبع + جع < وق 
)۲ 7 < 17 
)۳( ۲ - 32۶ 


در روابط بالا مه ضریب انبساط حرارتی حجمی و 
16 مدول بالک است و از رابطه‌ی (۶) به‌دست می‌آید: 


(( 7 ۷080 كِ 


پاش 
کون ان ۵ حجم مخصوص اولیه و 3 و ظ 
پاراترهایمدل: ی ناشند: 
هم‌چنین برای محاسبه‌ی کرنش برشی فقط از 
مژلفه‌ی کرنش پرشی مکانیکی استفاده شده و از کرنش 


برشی حرارتی صرف‌نظر شده‌است. 


6 مج + جع < 1 
اعاه 

1 < «8ع < 1۶ 

)1( چم و 2 


که در روابط (۵) و (1) و (۷) ظ پارامتر مدل» ۲ 
مدول بر شی در دمای اولیه و ۴ مدول بر شی در دمای 
ثانویه می‌باشد. 

در قسمت پلاستیک نیز باتوجه به اين که پارامتر 
سب خت‌شسو ندگی تنش پیش تحکیمی درنظر گرف ته 
ی ی سای ان ی( 


محاسبه شده‌است: 


۳ 1 
۸ تعك- - لول 
)۸( ۱ 
)۹( 1 
ی (0۵۰7- ود 


نشريه مهندسی عمران فردوسی 


مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


13 ۰( 

در روابط )٩(‏ و (۱۰) ۷۲ شیب خط حالت بحرانی 

و ۷ حجم مخصوص می‌باشد. هم‌چنین در این مدل 

رفتاری رابطه‌ای برای تبدیل 0 تنش پیش‌تحکیمی در 

دمای محیط به ۳,۲ تنش پیش‌تحکیمی در د مای بالا 
ارائه شده‌است. 


۲۲ < ۲۳ 90) 30۳ ۲0۵۲( (۱۱2 


پانفرنظر کف قاترن تهریان سای باب فادو از 
رابطه‌ی (۱۲ کرنش برشی پلاستیک محاسبه شله‌اسست. 


 _ 2 
ول‎ - 1 


_ ۵۳/2 )6-1(-۶[ 
۷ (۱۳ 


معادله‌ی سطح تسلیم نیز با اضافه شدن پارامترهای 


0 و 6 و 0 به معادله‌ی سسطح تسلیم مدل کم‌کلی 


اصلاح شده به‌صورت رابطه‌ی (۱ بازنویسی شده‌است: 


ِ ۳ -6) 02 
(و زر - ترچ ی م۲ - و 


20-1 


62 < ۱۵( 


برای نوشتن سابروتین یومت برای مدل رفتاری 
ملکور::ابعدا ماتر یس الاستیک و الانشتم پلاستتیک آن 
بهدست آمده‌انست که براین اباس برای تشکیل ماتریین 
ا لاس تیگ رارنظه‌ی( ۱ وی ای سکیا اما نی 
الامتوپلاس تیک مدل رفتاری از رابطه‌ی (۱۷) استفاده 


شده‌است [29]. 

3۱ 67 .(6). فك بت ۵ + 200 [*0] 
۳ (02(62)*« 

۱۷ ,19 مت 0۱ - ۱۱3 


۰۰8۵ (-() (ب) 


پس از تشکیل ماتریس‌هاء تابع تسلیم تعریف شده 
و شرط استفاده از هرکدام از ماتریس‌ها طبق تابع تسلیم 


برقرار شده‌است. 


سال سی و سوم شمارة سهء ۱۳۹۹ 


مهدی نحداپرست- علی جمالو - سیدعلی قریشیان امیری ۵۳ 
مطالعه‌ی پارامتری مدل رفتاری رابطه‌ی تنش تسلیم با دما؛ زسبت پیش تحکیمی و 


در این قسمت به‌منظور بررسی بیش‌تر مدل رفتاری 
تورچی [28] نمودار تتش کرنش به‌دست‌آمده با تغییر 
برخی از پارامترهای موجود در آن ارائه شده‌است که 
برنامه‌ی این مطالعه‌ی پارامتری در جدول (۱) آمده‌است. 


کاهش هریک از آنها به کاهش مقاومت برشی منجر 


شده‌است (شکل 6 


جدول ۱ برنامه‌ی مطالعه‌ی پارامتری مدل رفتاری تورچجی [28] بااستفاده از مدل عددی آزمایش سه‌محوری 


مطالعه‌ی پارامتری 
قسمت الاستیک تابع تسلیم دما 
‌ شماره‌ی مدل 
(طبع) | مه ا ظ | 6 | (وط0عظ | هل | ظ 0 ۸ | ۲ 
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۳۳0 ۳۳-200 2600(سسه 


شکل ۱ اثر تنش پیش تحکیمی و شیب خط حالت بحرانی بر نمودار تنش- کرنش خاک در دمای ۶۰ درجه: الف) شیب خط حالت بحرانی. 


ب) تنش پیش تحکیم 


سال سی و سوم. شمارة سه ۱۳۹۹ 


نشریه مهندسی عمران فردوسی 


و مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


تأثیر دما و پارامتر های موجود در رابطه‌ی تابع شده‌ا ست و بر خلاف آنها برای افزایش مقاومت بر شی 
تسلیم بر نمودار تتش- کرنش نیز در شکل (۲) مشاهده کاهش 0 مورد نیاز است. از سه پارامتر اضافه‌شده به 
باتوجه به شکل (۲) م شاهده می شود که افزایش بوده‌است. 
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ننش انحرافی(کیاو پاسکال) 
8 


100 
)2( 
0 
25 20 15 10 5 0 
کرتش محوری(درصد) 
۱2-4ع]سه. .1-جاع]سه 2-1)ع]سه 
500 


ننش انحر افی(کیاوپاسکال) 


)0( 


30 25 20 15 10 5 0 
کرتش محوری(درصد) 
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شتکل ۲ ار پارامتر 9و دما بر نمودان کتشی- کرنشن شاکه الف) پارامتر ۵ب دما 


نشریهُ مهندسی عمران فردوسی سال سی و سوم» شمارة سه. ۱۳۹۹ 


مهدی خحداپرست- علی جمالو - سیدعلی قریشیان امیری 


۵۵ 


۶ 8 


تنش الحر اف (کیلو پاسکال) 
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شکل ۳ اثر پارامترهای تابع تسلیم بر نمودار تنش- کرنش خاک در دمای ۶۰ درجه: الف) پارامتر 0 ب)پارامتر ۵ 


جدول ۲ مشخصات نمونه خاک براساس مدل رفتاری تورچی [28] در مدل آزمایش سه‌محوری 


71 2 ۵ 0 ۰ 


۸۳۴800 


85 025 0 014 3 18 


۵ ,2 1۷ داده‌های آزمایشگاهی 
۱ 

8 | ورن | ووج | چکرواک و لالوبی 
71 


# پارامتر 1 شیب خحط حالت بحرانی پارامتر 2 شیب خط تحکیم عادی. پارامتر ۸ شیب خط باربرداری و پنج پارامتر ۰7 ۰0 ۰1 » و 8 


پارامترهای مدل هستند. 


مدل‌سازی عددی آزمایش سه‌محوری (1) 
در مرحله‌ی اول برای بررسی عملکرد سابروتین یومت 
در مدلی با یک المان» مدل آزمایش سه‌محوری که 
هندسه‌ی نمونه‌ی مورد آزمایش در آن یک استوانه به‌قطر 
۵ سانتی‌متر و ارتفاع ۱ سانتی‌متر می‌باشد. ساخته 


شده‌است. 


سال سی و سوم» شمارة سه ۱۳۹۹ 


برای شبیه‌سازی رفتار خاک از سابروتین یومت 
استفاده شده و مشخصات خاک نمونه باتوجه به کد 
ارائه‌شده در حدول (۲) افله‌استت: 
شده‌است که در گام اول تحکیم خاک با تنش اولیه 


نشریه مهندسی عمران فردوسی 


۶ 


صورت گرفته و در گام دوم برای شبیه‌سازی فرآیند 
زهکشی, بارگذاری به‌صورت کنترل جابه‌جایی با سرعت 
تقریبی ۰,۶ میلی‌متر بر ساعت صورت گرفته‌است. برای 
مش‌بندی خاک نیز از المان هشت گرهی خطی (0318۳) 
که از گروه فشار آب حفره‌ای- تنش می‌باشد. استفاده 
شده‌است. 

تنش پیش‌تحکیمی خاک مورد بررسی ۸۰۰ 
کیلوپاسکال است که با فرض نسبت پیش‌تحکیمی ۲ و ۳ 
, تنش اولیه‌ی خاک به‌ترتیب برابر ۳۰۰ و ۲۰۰ 
کیلوپاسکال بوده‌است. هم‌چنین دمای اولیه ۲۲ درجه و 
دمای انویه ٩۰‏ درجه‌ی سانتی گراد فرض شده‌است. 


مدل‌سازی عددی شمع حرارتی 

برای شبیه‌سازی رفتار شمع حرارتی. هندسه‌ی شمع 
به‌صورت یک استوانه به قطر یک متر و ارتفاع ۲۹ متر 
(شکل ۶) و هندسه‌ی خاک اطراف آن به‌صورت یک 
توده‌ی خاک چند لایه به طول و عرض ۲۰ متر و ارتفاع 
۲ متر درنظر گرفته شده‌است. ضخامت و نام گذاری 
لایه‌های خاک مطابق شکل (۵) می‌باشد. 

مشخصات مکانیکی» حرارتی و هیدرولیکی 
لایه‌های مختلف خاک و هم‌چنین مشخصات مکانیکی و 


حرارتی بتن در جدول (۳) ارائه شده‌است. 


نشریه مهندسی عمران فردوسی 


مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 
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برای معرفی رفتار مکانیکی خاک لایه‌ی 9 علاوه‌بر 
مدل رفتاری دراکر- پراگر از مدل رفتاری تورچی [28] 
نیز استفاده شده که مشخصات آن در جدول (۶) 
آمده‌است. این مشخصات باتوجه به وجود نتایج صحرایی 
بااستفاده از روش کالیبراسیون حاصل شده‌است. 
پارامترهای جدول (۶) نیز مانند پارامترهای جدول (۲) 
تعریف می‌شوند. هم‌چنین برای معرفی مدل رفتاری 
کم‌کلی اصلاح‌شده تغییرات در سابروتین یومت مدل 
رفتاری تورچی [28] ایجاد شده‌است. باتوجه به شرایط 
اشباع خاک برای تحلیل شمع حرارتی از دو تحلیل مجزا 
استفاده شده‌است. در تحلیل اول برای تعریف آندرکنش 
حرارتی بین شمع و خاک از دو نوع اندرکنش ۲۳6۳۳0۵1 
۵ و 0606120108 ۲1621 استفاده شده‌است و 
در تحلیل دوم برای تماس مکانیکی بین شمع و خاک نیز 
از یک رفتار تانژانتی از نوع پنالتی با ضریب اصطکاک 
۰,۷ استفاده شده‌است. بارگذاری به سه صورت ازمون ۱ 
و ۲ و ۳ انجام‌شده که در آزمون ۱ فقط بار مکانیکی 
به‌اندازه‌ی ۱۳۰۰ کیلونیوتن به شمع وارد شده‌است. 
آزمون ۲ بارگذاری مکانیکی ندارد و فقط بار حرارتی 
به‌میزان ۱۳,1 درجه‌ی سانتی‌گراد و فرآیند تحکیم خاک 
درنظر گرفته می‌شود و درنهایت در آزمون ۳ بار مکانیکی 
و حرارتی و فرآیند تحکیم به‌صورت هم‌زمان اعمال 


سال سی و سوم شمارة سهء ۱۳۹۹ 


مهدی حداپرست- علی جمالو < سیدعلی قریشیان امیری 


شده‌است. دمای اولیه‌ی خاک نیز ۱۳,۱ درجه درنظر 
گرفته شده‌است. برای مش‌بندی شمع از المان‌های 
هشت‌گرهی خطی از گروه تتش سه‌بعدی با فن انتگرال 


105 
10 
10 
10 


10 


سال سی و سوم. شمارة سه ۱۳۹۹ 


جدول ۳ مشخصات شمع و لایه‌های خاک براساس مدل رفتاری دراکر- پراگر در مدل شمع حرارتی [18] 
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جدول ۶ مشخصات لایه‌ی سوم خاک براساس مدل رفتاری تورچی [28] در مدل شمع حرارتی 
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کاهش‌يافته (3181)) استفاده شده و المان‌های خاک نیز 
از المان‌های هشت‌گرهی خطی از گروه فشار آب 
حفره‌ای» تنش(318۳)) تشکیل شده‌است (شکل ۰). 


لایه‌های 


۵۸ 


+0006 ک0+ع5:00- 
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-3 0+06 
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مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


حش کتم(باکال) 


-3 506706 -3006+06 -2 50۴+ -2 06 


شکل ۷ نمودار مربوط به آنالیز حساسیت تأثیر ابعاد مش برای سه حالت مختلف (با تعداد المان‌های 6۵۰۰۰ ۵۰۰۰۰ و ۵۵۰۰۰) بر مقدار 


هم‌چنین در شکل (۷) نتایج آنالیز حساسیت تأثیر 
ابعاد مش. برای سه حالت مختلف (با تعداد المان‌های 
۰ ۵۰۰۰۰ و ۵۵۰۰۰) بر مقدار تتش قائم آمده‌است 
(تعداد المان‌ها رند شتسد ات : همان‌طور که مشاهده 
می شود تغییر سایز مش‌بندی تغییری در م سئله و روند 
آن ایجاد نمی‌کند. البته بدیهی است تغییرات قابل توجه 
در ا ندازه‌ی مش ها می‌تواند باعث عدم هم‌گرایی و 
ناپایداری در حل مسئله گردد. 
به‌منظور صحت سنجی مدل سازی‌های انجام شده نتایج 


نشریه مهندسی عمران فردوسی 


به‌د ست‌آمده از آنها با داده‌های صحرایی و آزمایشگاهی 
موجود در تحقیقات گذشته مقایسه شده‌است که 
جزئیات هر کدام قد فز سم آووفع مسق 
آزمایش سه‌محوری ۷ 

برای صحت‌سنجی سابروتین یومت نوشته‌شده‌ی 
نمودارهای تنش- کرنش و کرتش حجمی- کرنش 
محوری به‌دست‌آمده از مدل آزمايش سه‌محوری ) در 
شرایط ۲۲ درجه و ٩۰‏ درجه با نسبت پیش تحکیمی ۲ و 
۳ با داده‌های آزمایشگاهی گزارش‌شده توسط چکرواک 
و لالویی [7] مقایسه شده‌است که در ادامه نتایج مربوط 


سال سی و سوم شمارة سهء ۱۳۹۹ 


مهدی خحداپرست- علی جمالو - سیدعلی قریشیان امیری 


به نمودارهای تنش- کرنش ارائه شده‌است. در شکل 
(۸) نتایج با نسست پیش‌تحکیمی ۳ در دمای ۲۲ و ٩۰‏ 
درجه نشان داده شده‌است. 

مقایسه‌ی نتایج آزمایشگاهی و مدل عددی حاضر 
نشان می‌دهد که مدل رفتاری استفاده‌شده در تحقیق 
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کرش محوری(درصد) 


۵۹ 


حاضر توانسته در دمای ۲۲ درجه و ٩۰‏ درجه میزان 


دمای ۳۳ درحه قن فستهت ببستشوخ شون کی نمودار 


ضعف‌هایی مشاهده می‌شود. 


تنش اقحراقی «کیلوباسکال) 


نتایج ثحلیق حاضر سح داده های آزمایشگاهی(چکرواک و لالوبی [۷]) - -« 


15 20 25 30 


تتش انحرافی(کیلو پاسکال) 


کرنش محوری(درصد) 


نتایج تحقیق حاضر سس داده های آزمابشگاهی(چکرواکو لالوبی [۷]) > س 


شکل ۸ نمودار تنش- کرنش مدل عددی حاضر برای خاک با نسبت پیش تحکیمی ۳ ۵) دمای ٩۰‏ درجه. ) دمای ۲۲ درجه 


سال سی و سوم شمارة سه ۱۳۹۹ 


نشریه مهندسی عمران فردوسی 


مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 
نمودار تنش- کرنش خاک در دمای ۲ و ٩۰‏ حجمی- کرنش محوری به‌منظور صحت‌سنجی مدل 
درجه با نسبت پیش‌تحکیمی ۲ در شکل )٩(‏ آمده‌است. عددی حاضر برای خاک با ذسبت پیش تحکیمی ۲ و در 
باتوجه به شکل )٩(‏ می‌توان دریافت که سابروتین دمای ٩۰‏ درجه‌ی سانتی‌گراد در شکل (۱۰) آمده‌ا ست. 
نوشته‌شده دارای توانایی مناسب برای پیش‌بینی رفتار 


همان‌طور که مشاهده می‌شود در مجموع نتیجه‌ی 


به ابیت الم رضایت بخش می‌باشد. 


دمای متفاوت بوده اسست. هم چنین نمودار کرنش 


10 15 20 25 


کرتش محوری(ذرصد) 


نتایج تحالیق حاشر سس داده های آزمایشگاهی(چکرواک و لالویی [۷])- س 


تنش انحر افی (کیلوپاسکال) 
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کرنش محوری(درصدا 


نتایج تحقیق حاضر سس داده های آزمایشگاهی(چکرواک و لالوبی  *)]۷[‏ 


شکل ٩‏ نمودار تنش- کرنش مدل عددی حاضر برای خاک با نسبت پیش ‌تحکیمی ۲: ۵) دمای ٩۰‏ درجه «) دمای ۲۲ درجه 


نشريه مهندسی عمران فردوسی سال سی و سوم شماره سه, ۱۳۹۹ 


مهدی خحداپرست- علی جمالو - سیدعلی فریشیان امیری 


30 25 0 کرنش محوری(درصد) 
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2 ۲06۲۱۲۲6۲۱۲۵۱ سس 


۶۱ 
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کرنش حجمی (درصد) 


شکل ۱۰ نمودار کرنش حجمی- کرنش محوری به‌منظور صحت‌سنجی مدل عددی حاضر برای خاک با نسبت پیش‌تحکیمی ۲ و در دمای 


۰ درجه‌ی سانتی گراد 


صحت سنجی در تحلیل شمع حرارتی 
به‌منظور صحت سنجی عملکرد سابروتین نوشته شده در 
تحلیل شمع حرارتی, نتایج با داده‌های عددی و صحرایی 
گزارش‌شده توسط لالویی و همکاران [18] مقایسه 
شده‌است. 

در شکل (۱۱) مقادیر تنش قائم در طول شسمع 
تحت بار مکانیکی ۱۳۰۰ کیلونیوتن براس‌اس نتایج 
صحرایی و عددی گزارش‌شده توسط لالویی و همکاران 
[18] و هم‌چنین نتایج به‌دسست‌آمده از مدل حاضسر 
بااسستفاده از سه مدل رفتاری دراکر- پراگر» کم‌کلی 
اصلاح‌شده و تورجی [28] آورده شده‌است. 

همان‌طور که در شکل (۱۱) ملاحظه می‌شود 
هنگامی که شمع حرارتی تحت بار مکانیکی خالص قرار 
گرفته‌است. مقایسه‌ی نتایج عددی و صحرایی نشان 
می‌د هد که در لا یه‌ی 3 استفاده از مدل رفتاری 
ترمومکانیکی تورچی [28] و کم‌کلی اصلاح‌شده نتایج 
بهتری نسسبت به مدل رفتاری مکانیکی دراکر- پراگر 


سال سی و سوم. شمارهُ سه, ۱۳۹۹ 


داشته‌است. 

ق  )۱(‏ فا اناد عه دی طرن شم ده 
پا یان فاز گر مايش یک تحکیم ۲۸ روزه تحت بار 
بخرازقی ره رای تسطالعات رین و هب کاران 
[18] و نتایج مدل‌سازی حاضر نشان داده‌است. 

عطق فک (۱۲) که تن فاثم فسسمن #نفق با 
حرارتی را نشأن می‌د هد. می‌توان دریافت که مدل 
رفتاری تورچی [28] نسسبت به مدل رفتاری دراکر- 
پراگر و کم‌کلی اصلاح‌شده توانایی مناسب‌تری برای 
پیش‌بینی تنش قائم شمع دا شتها ست. هم‌چنین برخلاف 
نیمه‌ی اول لایه‌ی 18 در نیمه‌ی دوم آن مدل رفتاری 
تورچی [28] نتایج بهتری را حتی نسبت به مدل رفتاری 
کی کی که قرط الم و فاران 
[18] مورد استفاده قرار گرفته‌است ارائه کرده‌است. 

نتایج به‌دست‌آمده از مدل سازی حاضر و مطالعات 
لالویی و همکاران [18] مربوط به تنش قائم در طول 


نشریه مهندسی عمران فردوسی 


۶۲ مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


مکانیکی و ۱۳ در حه بار حرارتی در روز ۱۲ یک فرآیند تحکیم ۲۸ روزه در شکل (۱۳) آمده‌است. 
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داده های صحر ایی(لالوبی و هعکاران [۱۸])سه مدل کم کلی اصلاح شدوس 
نتابج عددی(لالویی و همکاران [۱۸])-0- 


شکل ۱۱ تنش قائم در طول شمع تحت بار مکانیکی ۱۳۰۰ کیلونیوتن 
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مدل دراکرپراگرسه هدل ترموهکانیکی در تحقیق حاضر -* 
داده های صحر آیی(لالوبی و همکاران [۱۸])-0- مدل گم کلی اصلاح شد :سب 
نعایج عددی(لالوبی و همکاران [۱۸])- 


شکل ۱۲ تنش قائم در طول شمع تحت بار حرارتی ۶ درجه 


نشریهُ مهندسی عمران فردوسی سال سی و سوم» شمارة سه. ۱۳۹۹ 


مهدی خحداپرست- علی جمالو - سیدعلی قریشیان امیری 


۶۳ 


تنش قائم شمع(مگاباسکال) 


2 3 4 


مدل دراکر پ راگرسو 


داده های صحرابی(لالوبی و همکاران [۱۸]) 6 


هدل ترموعکاتیکی در تحقبق حاضر سس 


مدل کم کلی اصلاح شدهسب 


نتابج عددی!(لالوبی و همکاران [۱۸]) " 


شکل ۱۳ تنش قائم در طول شمع تحت بار مکانیکی ۰ کیلونیوتن و بار حرارتی ۶ درحه 


به‌طور کلی بیشینه تتش قائم موجود در شمع یکی 
از مهم ترین پاراهرهای طراحی. شم‌ها هی‌با شلد بر این 
اساس پیش‌بینی بیشینه تنش قائم در تحلیل دقیق‌تر 
شمع‌های حرارتی دارای اهمیت ب سیاری است. باتوجه 
به شکل (۱۳) می‌توان دریافت که بيشینه تتش قائم شمع 
تحت بار گذاری هم‌زمان مکانیکی و حرارتی در لایه‌ی 
5 خحاک رخ داده‌است که استفاده از مدل رفتاری کم‌کلی 
اصلاح‌شده و تورچی [28] در اين لایه باعث افزایش 
قابل توجه هم‌خوانی نتایج عددی و صحرایی نسبت‌به 
تاش که لقع مراکم وراکن توافت رفتان ماک 
در این لایه را بر عهده دارد. شده‌است. 

باتوجه به قسمت اول صحت‌سنجی. می‌توان 
دریافت که مدل رفتاری ترموم کانیکی تورچی [28] 
به‌درستی به نرم‌افزار آباکوس اضافه شده‌است. نتایج 
موجود در قسمت دوم هم تأثیر قابل توجه مدل رفتاری 
خاک در بررسی رفتار شمع‌های حرارتی را تأیید می‌کند 
و ان سر فهد که اسظافه ان عدل رقارزه تیان کین 


سال سی و سوم. شمارهُ سه, ۱۳۹۹ 


تورچی [28] نسبت به مدل مکانیکی کم کلی اصلاح‌شده 
تغییری در نتایج به‌دستآمده ایجاد نکرده‌است. همچنین 
همانطرد که قیلا ااشانه ده آنالید ات هش اد 
مدل شمع حرارتی انجام شده و تغییر در اندازه‌ی مش‌ها 
تأثیر قابل توجهی در نتایج نداشته‌است. 


نتایج و بحث 

به‌منظور برر سی بٍ شتر تأثیر مدل رفتاری خاک بر رفتار 
شمع حرارتی و خاک مجاور آن. خروجی‌های مربوط به 
مولفه‌های قائم و شعاعی تنش و کرنش به‌وجودآمده در 
شمع و خاک مجاور آن در لایه‌ی ظ استخراج شده‌است. 

در شکل (۱4) تنش‌های قائم و شعاعی 
به‌وجودآمده در خاک مجاور شمع در پایان فاز گرمایش 
بار گذاری ترمومکانیکی آورده شده‌اسست. هم چنین 
کرنش‌های قائم و شعاعی به‌وجودآمده در خاک مجاور 
شمع نیز در شکل (۱۵) قابل‌مشاهده است. 

همان‌طرر کافی کل (۱۶ اه میا سیخ یی 


نشریه مهندسی عمران فردوسی 


۶۴ مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


مدل رفتاری خاک باعث تغییر تنش شعاعی و قائم خاک 
مجاور شمع حرارتی شده به‌نحوی که این تغییر در 
بعضی نقاط به ۱۰۰ کیلوپاسکال هم رسیده‌است. اما 
تفاوتی میان مدل مکانیکی کم‌کلی اصلاح‌شده و 
ترمومکانیکی تورچی [28] وجود ندارد. هم‌چنین باتوجه 
به شکل (۱۵) می‌توان تغییر میزان کرنش شعاعی و قائم 
خاک ناشی از تغییر مدل رفتاری خاک مجاور شسمع 
حرارتی را ملاحظه کرد. 

شکل (۱3۱) تغییرات تنش شعاعی در طول شمع 
حرارتی تحت بارگذاری هم‌زمان حرارتی به‌میزان ۱۳,۶ 
در حه و مکانیکی به‌اندازه‌ی ۱۳۰۰ کبلونیوتن را نشان 
داده‌است. 


مشاهده می‌شود همان‌طور که میزان تنش قائم 


لش شعاعی خاک"کبلوپاسکال! 
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مدل د. اک رکه 


مدل لرمومکلیکی 4 20 


لس مدل گم کلی اصلاح شده ۷ 21 
گ 


به‌وجودا مده در شسمع حرارتی برای خاک لایه‌ی 33 
بااستفاده از مدل رفتاری تورچجی [28] نسبت‌به حالت 
استفاده از مدل رفتاری دراکر- پراگر بیشتر شده در 
در همان لایه بااستفاده از کد یومت بیشتر از مدل رفتاری 
دراکر- پراگر به‌دست آمده‌است. هم‌چنین نتایج مدل 
مکانیکی کم‌کلی اصلاح شده تغییر قابل توجهی ذسبت به 

هم چنین میزان تأثیر مدل رفتاری خاک بر 
کرنش‌های شعاعی و قائم شمع نیز در شکل (۱۷) آمده 
که باتوجه به آن می‌توان دریافت که مدل رفتاری خاک 
مجاور شمع تأثیر ذسبتاً کمی در مولفه‌های کرنش شمع 


ا ر ‏ حست 
مدنی فرا گرب گر 


مدل کم کلی اصلاع شدو سوت 20- 


مدل ترمومکلیتی سن 21 


شکل ۱۶ تنش‌های قائم و شعاعی در خاک مجاور شمع در پایان فاز گرمایش بارگذاری ترمومکانیکی 
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شکل ۱۵ کرنش‌های قائم و شعاعی در خاک مجاور شمع در پایان فاز گرمایش بارگذاری ترمومکانیکی 


نشريه مهندسی عمران فردوسی 


سال سی و سوم شمارة سهء ۱۳۹۹ 


مهدی خحداپرست- علی جمالو - سیدعلی قریشیان امیری 
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شکل ۱ تنش شعاعی در طول شمع حرارتی در پایان فاز گرمایش بارگذاری ترمومکانیکی 
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شکل ۱۷ کرنش‌های قائم و شعاعی در طول شمع حرارتی در پایان فاز گرمایش بارگذاری ترمومکانیکی 


تحقیق عددی حاضر شامل مطالعه‌ی پارامتری روی مدل 
رفتاری ترمومکانیکی تورچی [28] و هم‌چنین بررسی 
عملکرد شمع حرارتی تحت بارگذاری حرارتی و 
مکانیکی بااستفاده از مدل‌های رفتاری مکانیکی و 
ترمومکالیکی»سی‌باشد. برای اشاله رفن سل .قاری 
ترمومکانیکی از قابلیت سابرویتن‌نویسی در نرم‌افزار المان 
محدود آباکوس استفاده شد و دو مدل سه‌بعدی ان 


سه‌محوری و شمع حرارتی تهیه شد. بر این اساس نتایج 


سال سی و سوم شمارة سه ۱۳۹۹ 


زیر تهفسنگا امله ات 

۱. در مدل رفتاری ترمومکانیکی تورچی [28] پارامترهای 
الاستیک تاثیر کمی بر تتشن تسلیم داردو پارامترهای 
اضافه‌شده به تابع تسلیم مدل رفتاری کم‌کلی 
اصلاح‌شده در رفتار خاک بسیار تأثیرگذار بوده‌است. 

۲ مدل رفتاری انتخاب‌شده برای خاک مجاور شمع 
حرارتی در تحلیل عملکرد شمع حرارتی مثر بوده و 
نتایج عددی حاصل از مدل رفتاری ترمومکانیکی 
تورچی [28] و کم‌کلی اصلاح‌شده نسبت به مدل 


نشریه مهندسی عمران فردوسی 


۶۶ مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


۶ مولفه‌های شعاعی و قائم تتش و کرنش خاک اطراف 


نزدیک تر بوده‌است. 


. به‌علت مقدار کم تنش به‌وجودآمده در خاک اطراف 


شمم و تغییر دمای کم خاک تفاوت قابل‌توجهی 
میان نتایج مدل ترموم‌کانیکی تورچجی [28] و مدل 
مکانیکی کم‌کلی اصلاح شده وجود ندارد. 


شمع حرارتی و هم‌چنین تنش قائم و شعاعی شمح 
حرارتی متأثر از مدل رفتاری خاک اطراف شسمع 
است. اما تغییر در مدل رفتاری خاک اطراف شمع؛ 
تغییر قابل توجهی در کرنش‌های شعاعی و فائم شمع 
ایجاد نمی کند. 
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8,)2004(۰ 209-22 ,۲0 

٩01 ۹۵۲001] ۰‏ طا 1۵۳۲0۵1۵00۵5 ۱۱۵۷2۲60 26 رهان ۵۶ دمتاهم۳0 لهتصدما۵عن0 ربا ,ات0 .8 
1995(۰) 5۱۷۵068 ,276 

تلقطوعض ۵۶۲ عصنام۷۲۵ م۲تابطتاعومت مطا تم رممصهر۴ اممتصمطمم‌طصرمهصی‌مط ]۳ ریله له بل .۲ مبصقصطعتیکا ...9 
2006(۰) ,36-45 ۲0۰ ب(6)1 ۷۵ رک6۵ 660۳۱۵ 0۴ تنل ۱۱۵۱۱۵0۵ و ها۲۵عطه) 

و 1656 طمتاه0زمعصمن ما 0۷۵۵وجا0) ق2 عصتاهع۲۱ رها 1۳06۵0 وردان ۵۲ معصق ماه ۰۲۷ وله اه ربا بقلق/۰10۳ .10 
1993(۰) ,170-1863 .۱0 ,(33)4 ,۷۵۱ ,۵۳۱0۵9۱۵۵106 دای 

8 گر 01275 4ماه۹20 ۲۵۲ ۱۷۲۱۵081 آهمنصمطهعمممممط] وتانمتافوم ۸ بلح اه بت .1 مافصتامظ ,11 
(1996) ,145-169 ۴0۰ ,(41)1 ۷۵۱ ,«وهآ0ع6 

,1 ۹0۲ 0۶ 6۰۵۷10۲ظ )صملجمم۵ ماع زمره 1[ مومتام۱۷]۵0 ,.۲ رهم۲۱ 220 ر.ل نک ره ۷۷ و2 نا ,۷/298 .12 


(2016) ,040160541-13 ۴0۰ ,(142)8 .۷۵1 ,۱۸۵6۱۵ ۱۵۱۱۵۵۲۱۱۱8 0۴ ۲۱۵1 ول 
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مهدی خحداپرست- علی جمالو - سیدعلی قریشیان امیری 


۷ 01275 طد قامع]ظ متهمم‌صو1 ۵ مصتاام۷۵0 وان مه ۷۱۵001160 یله باه ورگ مصقطقهتن 
(2001) ,608-621 ۲0۰ و(38)3 ۷۵۱ با 0 066016611۱6۵ 

و۳۵۵ 20 معصق مصصتان ۷ ۲۱۵0۸۵۵0 «المصصمط ۱ ره ره09222ظ 220 و.1 .1 ,200ع7عظ .]۷ .۳ ,دوه -مامضر 
2007(۰) ,144-154 .۳0 ,(89)1 ۷۵۱۰ ,زوماهع و۵۱۵۵ :0127 >اماعجهظ 5011 ۵۶ متناووع:۲ مج ۷۷ ۳0۲۵ 
01276 200۵۵4 ۲۵۲ ۷۲۵08 عجرم لومتصمطمصه-مصصممطآ مه وت رلم۲۳222 20 بیظ مرلته]۲ 
2010(۰) ,445-459 .0 ,(4)4 ۷۵ ۲۱۱۵۱۱66 ۲601۵6۱۱۱6۵ ۵ ول ۲۱۱۱۵۵0۵ 

01276 20264 عم ولع0م]۱۷ لومتصمطام‌مص-مصصنمط]. علض مصصمگ عم یلع موبلا .ظ مع۲08 
-2952 ۲0۰ و(37)17 ۷۵۱۰ رقع60۱۱۵6۵0) 1۱۱ ۱۷۵۱/۱0۵6 هک اه ۷۱۵۱۱۵۵ ۵۲ ال ۲۱۱16۵0 
2013(۰) ,2971 

۵ 012۷5 2۵۲۵۵0 ۲۵۲ م۷۱۵0 لهمنصمطهممصمممممط ممعیتاو- ۲۳/۵0 ۵ وله ا۵ .۷ ۰ظ ,ق۲۱۵8 
2016(۰) ,1059-1080 ۲۵۰ ,(40)7 ۷۵۱۰ 6607۱66۱۵۱ رز و۱/۵۱۵0 م۵ ۳۷۱۵۱۵۵۵ ۲۵۳ ۲۲۵ ول 
۶ ۲62۷10۲ مطا ۶ کصمتمعناوع۱۲] آممتمصتا از صرح اماصمصصته ۳۲ ربا رافتلاته ۷ 20 ,.۷۲ بطاتالط رما بتتا۵10] 
رکع۲60۲۱۱۵۵۱۵) ۱ ۵6 امه ۷۵۱۱۵ ۵ ۱۱۱۵۱۱۵۵۵0۳۵۵ ۳۱۱۵۰۱ ۲ممصه‌طم :۲ ۲۱۵۵۲ 2 
2006(۰) ,763-781 ۲0۰ ,(3008 ,۱۷۵ 

۵ ,۰ 0127 صز عمازظ برعتمحظ ۵۶ عبامز27ط0ظ لهم‌تممطمعصه-مصصمط 1 تم -عجمصا ریلة اه ,۷۰ .۲ ,111 
2020(۰) ,1-12 ۳0۰ ,(0)0 ,۷۵۱ رکت 066016612 

مج مدوم عصتععظ ۵۶ لاو لقمتمصصتا. رب رطعافعروظ مه ری بتصهنک ,۷ 00202725 
85 ر060)6000109) ۵2 صا ‏ فلژمگ 0مصتقرع-مصز۴ صز وملزظ برعتممرظ ۶ احمصمملتعو ممتاعلتآمعومن) 
2016(۰) ,)تن ,0.3 ,62016610 ۳۱۵۲8 ۵0/۵۲۵۱۱۵۵0 5۲۱۱۱۵۳۱۱۵۱۲۵0۵۵ 1۱ 0 

تقصصعهطامعن ۵۶ مد متتامصصهندظ مه فلورلممه لومنمهطهممومصصمط 1 وب رواتهطهعلقت 220 ریگ بناع928 
2017(۰) ,04017076 ۵۰( ,(17)9 ۷۵۱۰ ری6۵ 060۱۱6 ]۵ 0۵ ۲۱/۵۳۰۱۵۵۵۵ ص52 صا مملزظ روم« 
عمط 1 ما 0ماععزهاناگ 5011 ۲۱6-2۲۵60 0۶ معصدل مان ۷ عصتاام۱۷۱۵0 ۶0۲ مجمجح۸ مه وله ام ریم ب1۷]2 
2017(۰) ,896-901 ۲۵۰ ,(34)6 ۵ ۷ ۱۵ص 660۱۵6 00۵9 0۷۵1۵9 

۱۵۵۵۷10۲ [2عتصمامع]۱۷ مطا گم ممزهعتاعع ۲ص لهمتتمصصتا لا ریک رفتل0601812 220 ربق رومعافهل بر بکنکتت00ام5 ]1 
سل 6۵0 060ظ ک و( کل 0۵0۵01 و۱۵۲ ۱۱ رونام از دوعص لمح وملزظ «عتعممط مامت 0۶ 
.(569-575,)2016 ۳0۰ .05۸ 

۴ 0000 2 0۴ 1۵500886 مها ۶ تاه لهمزتعصتاط رما بتتاملقا 0ص2 ر.۲ ۸ بقتتصا عاام ره مخصطهونا زر[ 
۵ 0۲۱۴۵۵۱۵۲۶ ,۱۵205 لهمنصهدهعصمصمط.. ۵ فممتاممتطاصوم پجممتلزنا مهم فملزظ . رتم« 
2016(۰) ,126-142 ۴0۰ ,72 ,۷۵۱ رعت 660166120 

,(62)6 ,۷۵1 0601۵61۱۳۵۵6 ,8ع۱ز «عتمص ۶ منامز27ط0ظ لهعتصهطامعصه-مصنمط ]۰ یله اه ربا .ظ ب۵ ۷۲۸۲۷( 
2012(۰) ,503-519 .۳0 

صذّ دما «عمصظ ۵۶ مر مطاهظ آهمتصهط‌مصصمصصمط]. ربا .1 بتافتاظ مضه رب رتعطممهه بیظ .06 بطهتاه‌عاخ 
2014(۰) ,399-412 ۴۵۰ و(963 ۷۵1۰ 06۵0/۵61۵ ۸۵۵ ,122 اتمتاوه۱ظ ۲112 


۵۵۲ 0۴ 2۷10۲ظ0ظ لهمنصهطمع-م۵۱ ۱[ درم -مصمصا .1 ردتتعنع۲ 220 رب ۸ رعصه! و.؟ ۷۰ ,هناعل 
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سال سی و سوم شمارة سهء ۱۳۹۹ ان بای و عفر ویر 


۶۸ مطالعه‌ی عددی عملکرد شمع‌های حرارتی بااستفاده از مدل‌های رفتاری... 


٩200, ۸6۱۵ 06601۵61066, ۷۵۱۰ 12)4(, ۵۰ 729-737, )2017(.‏ 1۱۱ جر ع1زظ 
28 تورجی. سعید. "ارائه یک مدل رفتاری ترمومکانیکی برای خاکهای رسی غیر اشباع " پایان نامه کارشناسی ارشد. دانشکده‌ی فنی و 
مهندسی . دانشگاه خوارزمی تهران» (۱۳۹۲) 
۰ 0۱0102۷ع) ۵۴ رازه۲ نصا تو۵0) بصری! تصحتطما یلم دوم . مصنام0م۱۷ قمع واتمتامهاظ ام ربج .5 ب24[مصت0هه .29 
(1394) 285-307 


نشريه مهندسی عمران فردوسی سال سی و سوم شماره سه, ۱۳۹۹ 


